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3 Abstract 

o 

We determine the limit of the highest achievable photoemitted electron brightness, due to heating 
r S^ just after emission into vacuum, applicable to both relativistic and nonrelativistic photoelectron 



production. This heating is due to poorly screened Coulomb interactions equivalent to disorder 



induced heating seen in ultracold neutral plasmas. Owing to the multi-scale nature of the initial 

■^- interactions of such systems, an N-body tree algorithm is used to compute the universal scaling 

of the disorder induced heating in fully contained bunches, and is shown to agree well with a 

simple model utilizing the tabulated correlated energy of one component plasmas. We also present 

CO 

y—{ simulations for beams undergoing Coulomb explosion at the photocathode, and demonstrate that 

both the temperature growth and subsequent cooling must be characterized by correlated effects, 

X 

as well as correlation-frozen dynamics. In either case, the induced temperature is found to be of 



several meV for typical photoinjector beam densities, a significant fraction of the intrinsic beam 
temperature of the coldest semiconductor photocathodes. 



Great progress has been made in the cooling of photoemitted electrons, in particular 
by those working with negative electron affinity semiconductor photocathodes and groups 
pursuing electron emission from laser cooled atoms. Several semiconductor photocathode 
groups have measured photoemission of equivalent temperatures well below kT = 25 meV, or 
the thermal energy at room temperature, either via the cryogenic cooling of the photocathode 
PQ, or via the maintenance of a pristine photoemissive conditions in ultra high vacuum, 
allowing the low effective mass of conduction electrons to produce an effect equivalent to 
Snell's law in which the velocity spread of electrons is drastically reduced at the vacuum 
interface [2j. Furthermore, using laser cooled atoms which are photoionized, there has been 
production of electrons at temperatures equivalent to ~ 1 meV [3]. 

This temperature of emission, along with the accelerating field at the source E acc , are 
the only two parameters of the fundamental brightness (transverse, normalized) limit of 
photoinjectors [I], such that: 

_ mc*fe E acc 

BnAD ~ 2nkT (1) 

max 

where / is the bunch repetition rate, many next-generation photoinjector and ultrafast 
electron diffraction laboratories, whose applications depend critically on beam brightness, 
plan to merge such ultra cold photoemission sources with photoguns to produce relativistic 
particle bunches. For various applications, such bunches may contain from 10 6 — > 10 9 
electrons, of densities in the range of 10 17 — > 10 20 m~ 3 . 

The vast majority of all particle accelerator beam dynamics simulations and analytic 
calculations for relativistic photoinjectors rely on the space charge (SC) approximation to 
calculate the interparticle interaction, including nonrelativistic dynamics at the photocath- 
ode. Within the SC approximation, the Debye length of the electrons is assumed to be much 
larger than the average interparticle separation. Analytically, this allows the reduction of 
the full 6N dimensional particle distribution function pw( r i, Vi, r 2, ••-, r n , v n ) to a single par- 
ticle form p 6 (r, v). In beam dynamics simulations, this allows for the use of macroparticles 
representing many electrons, from which the self-consistent calculation and evolution of the 
mean field corresponding macroparticle distribution is made. 

However, for the lowest temperatures attained in semiconductor photoemission kT ~ 5 
meV [2j , and the lowest of the above densities, the Debye length is already on the order of the 
interparticle separation. Thus the space charge approximation is no longer valid. For larger 



densities, we expect the interaction lengths of the interparticle force to be much greater than 
the Debye length. Thus we expect some intrinsic velocity spread from the stochastic, poorly 
screened Coulomb forces. Further reduction of the photocathode emission temerature at 
sufficiently high density would be blurred out by this stochastic heating. In this Letter, we 
calculate the amount of this stochastic heating. 

Though certain effects of Coulomb collisions are well studied in beam physics, as in the 
Boersch effect and intrabeam scattering, this high density, low temperature limit represents 
a new regime for beam calculations. Electron collisions in beams are traditionally handled 
via the diffusive Fokker-Planck method [5] or via averaging over two particle interactions [5J, 
which both involve the reduction of the phase space dimensionality, as in the SC approxi- 
mation. Thus, such methods are inadequate to describe the fully-coupled initial conditions 
of a near-zero temperature bunch. 

The heating associated with the relaxation of a random, near-zero temperature distri- 
bution of charges is known to the ultracold neutral plasma (UNP) community. In such 
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systems, a cold gas is laser ionized, and after a time on the order of the r = 2-kuj, 

1/9 

2tt (noe 2 /m e eo) , the ion temperature is seen to rise significantly due to near-neighbor in- 
teractions, to nearly 2>kT c /2 = e 2 /Aneoa, where a is the Wigner-Seitz radius, a = (^/Attuq) 1 ^ 3 . 
In the traditional plasma physics terminology, this effect is referred to as disorder induced 
heating (DIH), and the effect is seen for plasma coupling parameters V = e 2 / 'Aire^akT of 
order unity or larger [7] . If we here assume a photocathode emission temperature of 5 meV, 
the above beam densities correspond to V = 0.2 — > 2. This estimate suggests that for the 
next generation of colder and denser electron beams, DIH serves as a hard limit to the 
achievable beam brightness. 

Disorder induced heating has been experimentally observed for the ions in a neutral 
plasma [8]. In cold neutral plasmas, the electrons equilibrate much faster than the ionic 
ojp 1 , and then serve to screen ion-ion interaction. An expression for the amount of DIH in 
ions was first given for Yukawa systems (for which electron screening, but not electron-ion 
recombination, is considered) in [9]. The initial ionic distribution is uncorrelated, and can 
be defined as the zero energy state. As the ions relax, the ion distribution begins to develop 
order, and the resulting correlations have an associated binding energy. This correlation 
binding energy can be calculated from the pair correlation function g(r), and is only a 



function of T and the electron screening parameter k = Aoebye/o: 

Uc _ jj = 1 e 2 n f g(r, V, n)d 3 r 
NkT ~ 2 47re J r U 

This binding energy is balanced by an increasing ion temperature. Calculating the increase 

in temperature thus only requires knowledge of U for a given electron screening, which have 

been tabulated via molecular dynamics (MD) simulations [9J. 

It has been demonstrated that a trapped, charged plasma, such as an electron bunch, 

is equivalent to the one-component plasma (OCP) model, in which the charges exist in a 

uniform neutralizing background pi)]. The external containing potential allows the initial 

uncorrelated state to be viewed as zero total energy. As electrons relax, correlations develop, 

and the presence of a "Coulomb hole" in g(r) for r < a creates an effective correlation binding 

energy equal to that of a one component plasma. In the case of an OCP, there is no second 

species to provide additional screening, and thus U is only a function of I\ Owing to this 

simplification, MD data has been fit to a power series relation for U = aT + bT 1 ^ 3 + c, 

where the coefficients are given in [10], a = —0.90, b = 0.590673, c = —0.26569. We may 

now write down the expression for the the final coupling Tf, and thus the final equilibrium 

temperature, for a fully confined, initially uncorrelated (Oj = 0), uniform distribution of 

charges with initial coupling Tj, analogous to that presented in [5]: 

Uf = l(^-l)=aT f + bTf + c (3) 

It is important to note the power series relation is quoted only for T > 1, however we 
find that for T < 1, the above expression is well approximated by Tj = Tf, or the limit 
of no DIH, as we expect. Eq. [3] is plotted in Fig. [I] For an ideal bunch with zero initial 
temperature, Tj — > oo, Eq. [3] can be solved to give Tf m 2.23. 

Practical relativistic photoelectron sources contain a number of complicating factors. 
Inside the emitting material (either a cold gas or crystal), the electron coupling, and thus 
DIH, will be dramatically reduced by the presence of the ionic/nuclear potential, and thus 
here we only need to consider DIH developing in vacuum. The beam density during emission 
will vary in time and with position across the bunch. The temperature associated with DIH 
should be reached on the order of u^ 1 after electron emission into vacuum, which for typical 
beam densities is on the picosecond scale. On this scale, we can neglect relativistic effects, 
such as the I/7 2 damping of the interaction force [5J. The coupling given by Tf = 2.23 = 



T eq , for a given local beam density, can then be viewed as the fundamental photoemission 
temperature limit. We consider the effects of a changing local beam density, as well as the 
nonrelativistic effects of acceleration during the evolution of DIH below, and we find it to 
correspond well with the prediction of Eq. [3j Rewriting it in terms of practical units, the 
heating induced in a zero temperature bunch is given by kT[eV] = 1.04 x 10 _9 (rao[m~ 3 ]) 1 / 3 , 
which is only a function of interparticle separation. Along with Eq. [TJ this forms the 
fundamental phoemission brightness limit. 

To test the prediction of Eq. [3j we have chosen to use a tree-algorithm electrostatic 
simulation package, with nonrelativistic, spherical, uniform particle bunches. Traditional 
particle-in-cell integrators have an intrinsic length scale — the spatial grid on which local 
fields are calculated. However, tree algorithms lack a spatial grid, and thus are better 
suited to our initial conditions, where both close interactions and long range correlations are 
significant. Several plasma treecodes have been developed, see for instance [TTj. However, 
considering the simplicity of the problem, we elected to directly modify a code originally 
written for gravitational interaction [T2] . 

We consider a randomly distributed spherical bunch of 10 5 particles with uniform density, 
and vary only the density (10 17 —¥ 10 20 m~ 3 ) and initial Gaussian velocity distribution (kT 6 
{0.25, 20} meV or kT = 0) , where we use the definition kT = m e ((vf) — (xiVi) 2 /(x 2 )), where 
m e is the electron mass, and i is a Cartesian coordinate. There is no breaking of spherical 
symmetry in our simulations, thus equipartition will always hold. For a fully contained 
uniform bunch, (xiv j) ~ 0, and the temperature is just a measure of average electron kinetic 
energy, with no spatial dependence. Containment is done in simulation via the application of 
a radial external force equal and opposite to the SC field. The parameters of the interparticle 
force calculation are the Barnes-Hut opening angle 9, and the leap frog integrator time step 
dt. The opening angle is the force calculation accuracy parameter that determines whether 
an electron is treated individually or lumped together with other similarly distant particles 
[j~3] . No effect was seen from a force-softening distance p, such that the force between two 
particles is /12 oc (r^ + p)~ 2 , from p = p = a/10 6 up to p = a/100. For all simulations 
presented, we have demonstrated convergence using dt = r/120, 9 = 23°, N = 10 5 . 

For an initial distribution with zero temperature, the final temperature for several ap- 
plicable beam densities is plotted in Fig. [2] We note that for extremely dense beams, 
n > 10 19 m~ 3 , the DIH is comparable to the temperature of photoemission from the 
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FIG. 1: Final coupling Fj, vs initial coupling Tj, given by Eq. v3\ (blue line), with Tf = Tj for 
Ti — > (red, dotted), and Tf = 2.23 for Ti — > oo. Dots are simulation results of a fully contained, 
equipartitioned electron sphere, density uq = 10 20 m -3 , and initial temperatures between kT=0.25 
meV and 20meV. Squares are simulation results for multiple densities, where kT z = 0, and kT x = 
kT y , and where T is calculated from the average of the three directions. Error bars are an estimate 
of the uncertainty in final temperature determination due to the residual oscillations, as in Fig. [2} 



coldest semiconductor cathodes kT ~ 5 meV. Furthermore, we note the oscillation of the 
electron temperature at 2uj p , which has been measured experimentally in UNP systems [8]. 
These oscillations last far longer than those seen in UNPs, as the bunch remains uniform in 
simulation (rather than Gaussian), and there is thus no spatial temperature smearing. 

For bunches with nonzero initial temperature, the amount of additional heating can be 
calculated in the final coupling, plotted along with the prediction of Eq. [3] in Fig. [TJ All 
simulations here had a constant density n Q = 10 20 m -3 , with varied initial Maxwellian 
velocity distribution. The duration of the temperature oscillations visible in Fig. [2] is re- 
duced with increasing initial temperature, though the exact dependence was not extracted 




FIG. 2: Temperature vs. number of plasma periods from treecode simulation using 10 5 contained 
electrons, and no = 10 20 (top curve), 10 19 , 10 18 , 10 17 m -3 (bottom curve), with zero initial 
temperature, Tj — )• oo. The equilibrium value predicted by Eq. [31 Tj w 2.23 is shown (dotted 
lines). 



in this study. In practical systems, initial acceleration in a voltage gap V will have cooled 
any photocathode emission temperature along the acceleration direction by a factor kTi/eV, 
making fcTji ~ almost instantaneously in high gradient electron sources, whereas the trans- 
verse velocity spread will remain unaffected [5]. Simulations were performed for contained 
bunches of multiple densities with kT x = kT y ^ and kT z = 0. Equivalent dynamics to 
Fig. |2j as well as agreement with Eq. [3] were found if the average of the temperatures in all 
three directions is used to calculate Tj and Tf, as shown in Figjlj The timescale of equipar- 
tition here also depends on kTi and no, but for the densities considered above and initial 
temperatures comparable to DIH heating, it was seen that equipartition occurs between r 
and 3t. A bunch instantaneously emitted into a uniform accelerating field will not have the 



DIH evolution altered, a fact that was also verified in simulation. 

Real electron bunches just after emission are very infrequently fully transversely contained 
in the sense described above, and are not emitted instantaneously. The overall space charge 
force can double the beam radius on the timescale of w~l, and acceleration can significantly 
lengthen the pulse afterwards. Thus, we must consider the process of DIH in bunches 
with time dependent density reduction. A system of fully coupled charges will have an 
equilibrium temperature that decreases as kT oc 1/T eq a, as given by above. However, 
it is well known that system of noninteracting charges expanding under linear forces will 
have an equilibrium temperature that decreases as kT oc 1/a 2 [5J, as a consequence of 
RMS emittance conservation, which therefore leaves the phase space density or brightness 
unchanged. During the explosion, the beam is in a highly nonequilibrium state, but our 
definition of temperature above still applies, where here (rv r ) ^ 0. 

Simulations were performed to determine the temperature as a function of time for a 
bunch undergoing full Coulomb explosion. Multiple densities were considered, however the 
previous scaling with n was seen. Thus, we present only results for n = 10 20 m -3 in 
Fig. [3J The total overall expansion and velocity growth is plotted in the inset, with the 
analytic prediction, showing excellent agreement. The temperature shows an initial increase 
due to DIH at a time of £ max ~ u p 1 , to a value of kT(t = £ max ) = e 2 /47reoa(£ max ), or a value 
a(£ max )/a ~ 1.6 times smaller than for a fully contained beam. The bunch continues to 
expand, attempting to equilibrate to T eq = 2.23, and the temperature continues to decrease. 
However, at longer times, the temperature falls as a -2 , as in an uncoupled system. Thus, 
there must be a time at which correlation ceases to grow with decreasing temperature. Such 
behavior has been noted in the adiabatic expansion of UNPs due to kinetic pressure [H] . 

Thus, to model the temperature as a function of time, we can write the temperature as 
a product of the correlated growth, and the decoupled expansion after correlation ceases. 
Since the growth of the radius is small during the time of DIH, we may presume T = 

2 2 

T c {t)-^fy2, where T c is the correlated temperature dependence, and the term -^%z describes 
the expansion of independent particles. The rate of DIH should be proportional to the how 
far the system is from the equilibrium temperature at the current density, and so: 

*E° = T ^ ~ T ^ = 1 ( t T(t)] (4) 

dt r r V4vre a(t)r e(? cy ' J y ' 

The only free parameter of this model is decoupling time To, or the time when the corre- 

8 



lation ceases to increase. Integrating this numerically, and fitting to the temperature data, 
we find r = 0.325r, which is also plotted in Fig. [3} This time is non-negligibly greater 
than t max , and suggests that the bunch is cooled not only via decoupled expansion, but 
coupled expansion as well just after t max . We can verify this prediction by looking at the 
pair correlation as a function of time. This is computed in Fig. [4], via averaging over 5 x 10 3 
particles randomly selected from the distribution. The development of the "Coulomb hole" 
and an accompanying shock profile due to violent repulsion is clearly visible, as has also 
been seen in UNP simulations [15]. Indeed the correlation ceases to change after a time of 
approximately r = 0.3r. Partial confinement effected by linear focusing in the source would 
yield an altered a(t), and thus a longer To, whereas acceleration that lengthens the bunch 
significantly near the DIH heating timescale would yield a shorter tq, but in either case the 
DIH heating timescale should remain unchanged. 

The exploding beam case justifies post-factum the validity of the temperature limit of 
Eq. [3] for practical photoemission, though beams in many applications do not have constant 
volume. In the exploding case, near r the temperature scaling shifts from a(t) _1 as given in 
Eq. 2, to the familiar scaling of a 2 oc a(t)~ 2 , as given by emittance/brightness conservation. 
A more detailed calculation of DIH for practical cases in which the bunch shape is both 
time and space dependent should involve averaging over a temperature distribution given 
by Eqs. [3] and [4] using the local density as determined by space charge tracking. 

In summary, this Letter has proposed that the fundamental temperature limit of pho- 
toemission will be from the disorder induced heating of electrons due to poorly screened 
Coulomb interactions at all length scales. We have shown that the tabulated thermody- 
namic quantities of one component plasmas are sufficient to explain both fully contained 
and Coulomb exploding instances of DIH, and have verified two simple relations that de- 
scribe the temperature evolution. Furthermore, many interesting effects of DIH in UNP, 
such as temperature oscillation and correlation decoupling, should also be present in such 
cold beams, yielding the possibility of rich interdisciplinary study. 

Practically, for next-generation ultracold dense electron sources we have shown that the 
temperature of photoemission, and thus the maximum beam brightness, cannot be arbitrar- 
ily improved. Furthermore, given the rapid progress of photocathode temperature reduction, 
we anticipate such a limit to be approached in the next generation of high brightness electron 
sources producing intense beams. 
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FIG. 3: RMS velocity spread of a spherical bunch undergoing coulomb expansion, treecode 
data (circles), and prediction of Eq. [4l The inset shows the overall expansion of the bunch in 
normalized units, both analytic prediction (solid line), given by Gauss's Law, and treecode data 
(dots). Uniform bunches remain uniform, and the normalized bunch size R{t)/Rq is independent 
of density. 
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